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WICK PRODUCT IN THE STOCHASTIC BURGERS EQUATION: A 

CURSE OR A CURE? 

S. KALIGOTLA AND S. V. LOTOTSKY 

Abstract. It has been known for a while that a nonUnear equation driven by singu- 
lar noise must be interpreted in the re-normalized, or Wick, form. For the stochastic 
(<— ^ I Burgers equation, Wick nonlinearity forces the solution to be a generalized process no 

psj ■ matter how regular the random perturbation is, whence the curse. On the other hand, 

certain multiplicative random perturbations of the deterministic Burgers equation can 
? ' only be interpreted in the Wick form, whence the cure. The analysis is based on the 

study of the coefficients of the chaos expansion of the solution at different stochastic 
t^^ ■ scales. 

0_i ■ 1. Introduction 

'^ ! Burgers equation, 

(1.1) Ut + uUx = Uxxi t > 0, X G M, 

first suggested as a simplified model of turbulence (Bateman [T], Burgers [2] E]) is 
now used to study problems such as traffic flows (Chowdhury et al. [4j) and mass 
^€1 • distribution of the large scale structure of the universe (Molchanov et al. [5]). The 

Q\ '. equation also appears in the study of interacting particle systems (Sznitman [6]). 

^ I Considering random initial condition and/or driving force in equation (11. ip is one 

• ' way to model and investigate turbulence. The idea goes back to Burgers himself and 

Q ■ the result is known as the Burgers turbulence. A popular option for the driving force is 

O ! additive space-time white noise W{t, x). Mathematical theory of the resulting equation, 

(1.2) Ut + UUx = Uxx + W{t,x) 

K> ! has been developed (Bertini et al. [7]). A more general version of (II. 2p with multi- 

H ' plicative noise, 

(1.3) Ut + UUx = Uxx + fit, x) + {g{u{t, x))^ + h{u(t, x))W(t, x), 

has also been studied (Gyongy and Nualart [8j). It turns out that (11.20 cannot be 
generalized much further while staying within the same mathematical framework, when 
the solution of the equation is a square-integrable random variable with sufficiently 
regular sample paths. In particular, equation 

(1.4) Ut + UUx = Uxx+UxW{t,x) 

makes no sense, because the solution, if existed, cannot be regular enough to define 
the point-wise multiplications uUx and UxW. 

Using the Wick product o instead of the usual point-wise multiplication makes it 
possible to study (ll.4|) and similar equations in the framework of white noise theory. 
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The operation o was first introduced in quantum field theory (Wick |9]). In stochastic 
analysis, the operation is essentially a convolution (Hida and Ikeda [ID]), and is closely 
connected with the Ito and Skorokhod integrals (Holden et al. [TT], Chapter 2]). The 
Wick version of (11.41) . 

(1.5) Ut + UOUx = Uxx + u^o W{t, x), X EM., 

has been studied (Benth et al. [121 Section 5.2, Example 2]), and is known to have a 
solution in the space of generalized random processes. The tools of white noise analysis 
and availability of generalized random processes make it possible to consider even more 
singular noise than W, both additive (Holden et al. [131 E]) and multiplicative (Benth 
et al. [12]). 

Our objective in this paper is to study the equation 

Ut{t, X) + U{t, X) O Ur,{t, X) = Uxxit, x) + fit, x) 

(1-6) + ^ {ak{t, x)uxx + h{t, x)ux + Cfc(t, x)u + gk{t, x)) o ^k, 

k>l 

<t <T, u{0,x) = ip{x), where {^k, k > 1} are independent identically distributed 
(iid) standard normal random variables, T < oo is non-random, the coefficients ak,bk, Ck 
and the free terms /, Qk are non-random, and x E G, with 

• G = M (whole line) or 

• G = S^ (the circle, which corresponds to periodic boundary conditions). 

Equation (11.61) includes (11.51) as a particular case because the space-time white noise 
W(t,x) can be written as 

(1.7) W{t,x) = J2hk{t,x)^k, 

k>l 

where {hk, k > 1} is an orthonormal basis in L2((0,T) x R) (see, for example, Holden 
et al. [TTl Definition 2.3.9]). From the physical point of view, (II. 6p is a natural random 
perturbation of the original equation (11.11) . Indeed, one possible interpretation of (11.11) 
is the motion of a one-dimensional fiuid, and then the basic laws of fiuid dynamics 
suggest a more general version of (II. ip : 

(1.8) ut + uux = -{fJ^Ux)x + -F{t,x), 

p p 

where p is the density of the fiuid, n is the (dynamic) viscosity, and F is the external 
force (see, for example, Gallavotti [I5| Section 1.2.2]). Thus, (II. ip is (ll.Sp with constant 
/J, and p. If p is not known, then (assuming p is still constant) a possible approach is 
to consider 

(1.9) p{t,x) = po + eW(t,x), po = const; 

given the time and space scales of the model, one can choose e > small enough 
to have the right-hand side of (II. 9p positive with probability arbitrarily close to one. 
Substituting (II. 9p into (II. 8p and using (II. 7p (and interpreting all multiplications in 
the Wick sense) leads to a particular case of (II. 6p with ak = po + hk, hk = dhk/dx, 
f = Ck = gk = 0. Note that we must interpret multiplications in the Wick sense: no 
matter how small the e is in (II. 9p . the equation ut = {I + eW(t,x))uxx is ill-posed 
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path- wise. An extra benefit of this interpretation is preservation of the mean dynamic, 
something we would never get with the usual Burgers equation. Indeed, a remarkable 
property of the Wick product is that the (generalized) expectation of the product is the 
product of expectations: K{uov) = {Eu){Ev). As a result, the (generalized) expected 
value U = Ku of the solution of (II. 6p satisfies the usual Burgers equation 

Ut + UU, = U,, + /, U\t=o = Eyp. 

The main results of the paper are as follows: 

(1) If the functions ak,bk,Ck, f,gk are non-random, bounded and measurable and 
the initial condition </? is a generalized random field on G, then (under some 
addition technical conditions on / and ip) there exists a unique generalized 
process solution of (11.61) : 

(2) The Wick product uoux forces the solution of (I1.6P to be a generalized process 
even when the random perturbation is very benign (such as at = h^ = Ck = 
gk = for all /c > 1, and ip{x) = ^i h{x) for a smooth compactly supported h). 

Our approach to proving existence and uniqueness of solution is to derive the chaos 
expansion of the solution of (II. 6p and to establish explicit bounds on the coefficient of 
the expansion. The proof has a combinatorial component due to the appearance of the 
Catalan numbers. Mikulevicius and Rozovskii [16] use the same approach to study the 
Wick-stochastic version of the Navier-Stokes equations. 

The overall conclusion, which also explains the title of the paper, is that, while so- 
lutions of any stochastic Burgers equation in the Wick form are generalized random 
processes (even when the use of the point-wise multiplication leads to a classical so- 
lution — whence the curse), certain random perturbations can only be considered for 
the Wick form of the equation (because point-wise multiplication cannot be defined 
— whence the cure). In our opinion, the cures (the ability to consider more gen- 
eral stochastic perturbations, preservation of the mean dynamics, and an easy access 
to the chaos coefficients of the solution) outweigh the curses (necessity to work with 
generalized random elements and questions about physical interpretation of the model). 



Although there are some similarities between our work and that of Benth et al. 
there are also several important differences: 

(1) The random perturbations in (II. 6p and in |12j are different and there is no 
obvious way to reduce one to the other; 

(2) The solution of (II. 6p is global in time (can be constructed for all T > 0); 

(3) The solution of (II. 6p is constructed on an arbitrary probability space, not just 
on the white noise space; 

(4) Our analysis of the chaos expansion of the solution leads to more detailed 
information about the solution space and can be used to compute the solution 
numerically. 

In Section [2] we outline the main constructions in the theory of generalized processes. 
In Section [3] we define the chaos solution for equation (11.61) and state the main result 
about existence and uniqueness of solution (Theorem 13. 3p . We also show (Theorem 
13. 2p that the Wick product indeed forces the solution to be a generalized process. In 
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Section m we derive the propagator for equation (11 .Op (the system of deterministic par- 
tial differential equations describing the propagation of random perturbation in the 
equation through different stochastic scales) and outline the proof of Theorem 13.31 We 
also discuss briefly the general technical issues related to the study of (11.61) and similar 
equations. The details of the proof of Theorem 13.31 are in Section |5l 

Throughout the paper, we fix a probability space F = {Q, J^, P) and a countable 
collection of independent and identically distributed (iid) standard Gaussian random 
variables ^ = (^1,^2, • • •) on F. We also assume that J-" is generated by $,. By R and 
C we denote the sets of real and complex numbers, respectively; C and C^ denote the 
spaces of bounded continuous and bounded continuously differentiable functions (with 
the derivative also bounded). 

2. Gaussian polynomial chaos and the Wick product 

We start with a review of some constructions from the white noise theory. 

Definition 2.1. A generalized Gaussian chaos space is collection of the following four 
objects (¥,^,H,q): 

• A probability space F = {Q,J^,F); 

• a collection of iid standard Gaussian random variables ^ = (^1, ,^25 • • •) such that 
T is generated by ^; 

• A separable Hilbert space H ; 

• An unbounded, self-adjoint positive- definite operator Q on H such that Q has 
a pure point spectrum: 

Qf)fc = qki)k, k>l, 

the eigenf unctions f)^ o/Q form an orthonormal basis in H , and the eigenvalues 
qk ofQ satisfy 

(2.1) y:4<oo 



fc>l ^k 



for some 7 > 0. 



Condition 
nuclear space 



2.ip ensures that the projective limit of the domains of Q", n > 1, is a 
To simplify some of the future computations, we will assume that 



(2.2) 1 < gi < g2 < • • • and V — < 1. 

There is no loss of generality involved, as we can always replace the original operator 
Q with cCf , with c sufficiently large. By q we denote the sequence (gi, g2, • • •) of the 
eigenvalues of Q. 



In the white noise approach, the dual of that nuclear space with the Borel sigma-algebra and a 
Gaussian measure is the probability space. We are able to use an arbitrary probability space without 
any topological structure. For us, the numbers qk are just weights to ensure convergence of certain 
infinite sums. 
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Example 2.2. (1) The traditional white noise constructions correspond to H = L2[ 
and Q = —A + |xp + 1. If rf = 1, then q^ = 2k. (2) Stochastic partial differential 
equations in R'^ correspond H = L2((0,T) x M"'), Q = -d^/dt"^ - A + |xp + 1, with 
periodic boundary conditions in time. (3) Stochastic partial differential equations in 
domains or on a manifold correspond to H = L2((0,T) x G), where G is a smooth 
bounded domain or a smooth compact manifold in M'^, and Q = —d'^/dt^ — A, with 
periodic boundary conditions in time and appropriate boundary conditions for the 
Laplace operator A on G. 

Next, we review some of the notations related to multi- indices. 

A multi-index o: is a sequence ex. = [ai,a2, . . .) of non-negative integers, such that 
only finitely many of a^ are different from 0. The collection of all multi-indices is 
denoted by J'. By definition, 

• /3 < a if /Sfc < a/c for all k; 

• /3 < a if /3fc < Ofc for all k and f3k < dk for at least one k. 

If /3 < a, then a — /3 is the multi-index {ak — (3k, k > 1). For a. E J define 

\CX\ = ^Ofe, Oc\ = JJctfe!. 

k k 

Special multi-indices and the corresponding notations are (i) (0) , the multi- index with 
all zeros: |(0)| = 0; (ii) e(A;), the multi-index with 1 at position k and zeroes elsewhere: 
|e(fe)| = 1. 
For a sequence z = [zi, Z2, . . .) oi complex numbers and ex. E J^ we write 

z" = Ylz'^", z'" = {z"y, reR. 

k 

Here are some useful technical results. 

Proposition 2.3. Let q = (gi,g2,---) be the sequence of eigenvalues of the operator 
Q. Under the assumptions (12.21) . 

(2.3) 

(2.4) 

(2.5) 

Proof. Define pk = l/qk- To establish (12. 3p and (12. 4p . note that, since limfc-^ooPfc = 0, 
we have 

Ok 



^ a! 


= exp 


\k>i 


1 

Qk 




A:>1 ^' 


Qk 

c-1' 


) 


a|! <q° 


a\. 







exp 



k k aeJ k '^ k>l ^'^ aeJ k 

To establish (12. 5p . let n = \ex\ and use the multinomial formula and (12. 2 p to find 

V k ^"Z aej';|a|=n ^ 

This concludes the proof of Proposition \2.3\ . D 
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By a theorem of Cameron and Martin [TTj, every square- integrable function -F(^) 
with values in a (complex) Hilbert space V can be written as 

(2.6) F(^) = J2 F»U 

where 



H„,(x) = (-l)«e^V2^,-.V2 

and 

F« = E(F(|)e„). 
Then 

E||F(Olly=Ell^-llv- 

We denote by L(^; V) the collection of all square-integrable l^- valued functions -F(^). 
Next, we construct spaces of stochastic test functions and generalized random ele- 
ments. 

Definition 2.4. For p e [0, 1] and i>0, 

• the space {S)p^i{V) is the collection of ^ E IL2(^; V) such that 

aeJ 

• the space {S)^p^-e{V) is the closure o/L2(^; V) with respect to the norm 

ll^l|2 V^/'^ n-p.,-^c»iirT> l|2 

aeJ 

• the space {S)p{V) is the projective limit (intersection endowed with a special 
topology) of the spaces {S)p^£{V), as i varies over all non-negative integers; 

• the space {S)-p(y) is the inductive limit (union endowed with a special topology) 
of the spaces {S)-p^-e{V), as i varies over all non-negative integers. 

A stochastic test function is an element of {S)p{<C) for some p > 0. 
A (V -valued) generalized random element is an element of {S)^i{V). 

Thus, every $ G {S)-i{V) has a chaos expansion 
(2.7) <!> = J2^c.U 

a 

and the coefficients ^a provide information about $ at different stochastic scales (that 
is, in the linear spans of ^q, |q:| = n). 

In the white noise theory, the spaces (5)_o(C) and (iS)_i(C) are known, respectively, 
as the spaces of Hida and Kondratiev distributions (note that indeed (5)_o is not the 
same as (i5)o). 

It follows from 02.5|) that a\ in the definitions of the spaces can be replaced with |a|!: 
it will not change {S)p and (5)_p and will shift the index i in the individual {S)±p^±i. 
We will see below why the values of p are restricted to [0, 1]. 
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For ^ G {S)-p{V) and r] G (5)p(C), define (*,r/) G V by 
(2.8) (vl>,r/) = 5^vl/^r7«. 

By the Cauchy-Schwarz and the triangle inequahties, 

(2.9) ll(^,r])lk< (Eii*"IIv("0-V") ( E I^"I'(«')V" ) ■ 

Let z = {zi,Z2, . . .) be a sequence of complex numbers such that J2k>i \^k\'^ ^ °°- 
The stochastic exponential of z is the random variable 






Direct computations using (12. 3p . (12 ■4p and the generating function formula for the 
Hermite polynomials show that 

(1) If Efc^fckfcP < oo, then S{z) G (5)p,,(C) for every < p < 1; 

(2) Efc?fkfcP<l>then£(z)G(5)i,,(C). 

Definition 2.5. The S-transform ^ of ^ E {S)-p-e(y) with chaos expansion (12. 7p is 

$(^) = ($,£:(^)) = 5^ '^- ~- 






defined for z such that £{z) G (iS)p^£(C). 



Note that $(o) = $(0); $(o) is called a generalized expectation of $ and even denoted 
by E$, although it is clear from the corresponding definitions that a typical element 
of (iS)_p._£ for p,i > cannot have any moments in the usual sense. 

The following is a (partial) characterization of the spaces S-p{V) in terms of the 
S'-transform. 

Theorem 2.6 (Characterization Theorem), (a) // $ G (5)_p,_£ and < p < 1 then, 
for every real sequences p and r, such that £{p) G (5)p^^(C), £{r) G (iS)p^£(C), i/ie 
function f{z) = ^{zp + q) is an entire function of z E C 
(b) For < R, i <oo let 

(2.10) KeiR)=lz : ^ q^"|z"p < i?^ 

//$ G (iS)_i(V"), it/ien t/iere exist R,i such that ^{z) is analytic in Ke{R). Conversely, 
if f = /(^) ^s fl function analytic in Ki{R) for some < R,i < oo, then there exists 
a unique $ G (iS)_i(y) such that <!> = /. 

Proof, (a) See Kuo [TS^, Theorem 8.10]. There is also a converse statement character- 
izing entire functions with certain growth at infinity as ^-transforms of elements from 
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(b) See Holden et al. |11> Theorem 2.6.11]. Note that, given / = $, we recover 
(2-11) ^c = ^= ^..^1 



^\dz1'dz^'' 



z=0 

D 



In what follows, we refer to 'K^{R) from fl2.10p as a neighborhood of 2 = 0. 

The significance of Theorem 12.61 is that it provides an intrinsic characterization of 
the spaces {S)-p{V) in terms of the S'-transform and the spaces of analytic functions, 
as opposed to the extrinsic characterization of Definition 12.41 using a particular or- 
thonormal basis in L2(^; V). We also see why the values of p are restricted to [0, 1]: for 
p > 1, stochastic exponentials are never test functions, making it impossible to define 
the S-transform. 

Finally, we define the Wick product o for two elements of (iS)_i(C). 

Definition 2.7. Given $, \E' G (5)_i(C), $o\I^ is the unique element o/(5)_i(C) such 
that 

(2.12) $^(2) = $(z)$(;2) 

An immediate consequence of (I2.12p is a useful property of generalized expectations: 

E($o^) = (E<I>)(E^). 

Formula (12.121) is motivated by the following property of the Wick product (which can 
be traced back to the original paper by Wick [9]): 

(2.13) H„(6).H„to) =/""■«'='""'«''• "^J- 

Writing Hq,(^) = y/a\^ct, we conclude from (I2.13P that 

(2.14) H„(0oH^(0=H,+;3(0. 
If 

a ct 

then formal term-by-term multiplication using (I2.14p implies 

(2.15) $o^ = ^<l„^^H„+^(^), 

a,l3 

which is consistent with (I2.12p . Note also that if 
then \l/„ = <l>„\/al and 



(2.16) $ o M/ = 5^ ^l^_j^$^vl;^^^^^, 

which is not as convenient as (12.151) . 
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3. Existence and uniqueness of solution 

Let G = W or G = S^ (the circle, which corresponds to periodic boundary conditions 
and which we identify with R/(0,7r)). Denote by H'^{G), 7 G M, the Sobolev spaces 
on G and by || • ||^ the corresponding norms; H'^IG) = L2{G). Recall that, for G = R, 



'l + y'y\f{y)\'dy, 



where / is the Fourier transform of /, and for G = S^, 

k 

where fk are the Fourier coefficients of /. Also recall that, by the Sobolev embedding 
theorem, every element (equivalence class) from H^{G) has a representative that is 
a bounded and Holder continuous function on G. To avoid unnecessary technical 
complications, we will simply say that if / G H^{G), then / is bounded and Holder 
continuous, and sup^g^j 1/(3^)1 < C*!!/!!! for some G independent of /. 
Consider the equation 

Ut{t, x) + u(t, x) o u^(t, x) = u^x{t, x) + f{t, x) 

(3.1) + ^(afc(i^, x)u^^ + hit, x)u^ + Cfc(t, x)u + gk{t, x)) o ^k, 

k>l 

<t<T, X eG, 

with initial condition m(0, x) G {S) ^i{L2{G)) , where C,k are iid standard normal random 
variables, and /, a^, bk, Ck and Qk are non-random measurable functions. To define the 
solution we assume that f,gk G L2((0,T) x G) and 

(3.2) sup|afc(t,a;)| + sup\bkit,x)\ + sup |cfc(t,x)| + ||fi'fc|U2((o,T)xG) < ql 

for Qk from Definition 12. II and some r > 0. 

Definition 3.1. The process u G (5)_i((L2((0,T); if^(G))), is called a strong chaos 
solution of (13.11) if there exist < R,i < 00 such that, for all z = {zi,Z2, . ■ ■) G ^e{R) 
and all t G [0,T], the equality 

u{t,x;z)+ / u{s,x; z)ux{s,x; z)ds = u{0,x; z) + / Uxx{s,x; z)ds 

(3.3) ^ , . 

+ Z_. [ak{s,x)uxx + bk{s,x)ux + Ck{t,x)u + gk{t,x)]zkds 
k>i-^o ^ ^ 

holds in L2{G). 

The a priori assumption u G (5)_i((L2((0, T); //^(G))) implies 

{u,r])eL2{iO,Ty,H\G)) 

for all rj G (iS)i(C). The 5'-transform -u of m is thus an element of L2((0,T); H'^{G)), 
so that Ux G L2{{0,T);H\G)) and Ux e L2((0,T) x G). By the Sobolev embedding 
theorem, u,Ux G L2((0,T); C(G)), making the product uUx well defined. Condition 
(13. 2p ensures uniform in (t, x) convergence of all infinite sums. 
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The chaos solution is a variational solution in the space (iS)_i(L2((0, T); if^(G'))), 
and the characterization theorem (Theorem 12. 6p makes it possible to restrict the set 
of test functions to stochastic exponentials. The resulting deterministic equation for 
the 5'-transform ( 13. 3p (which follows from ( 13. ip by the linearity of the S'-transform 
and the definition of the Wick product (12.120 ) can be satisfied in a variety of ways 
(classical, variational, viscosity, etc.) Our a priori assumptions on u allow us to satisfy 
the equation for u in the strong sense, because all the necessary partial derivatives of 
u in t and x exist as L2 functions. This is why the resulting solution of (13. ip is called 
a strong chaos solution. 

While our main objective is to establish existence and uniqueness of solution of (13. ip . 
we will start by showing that the space (iS)_i is indeed the natural solution space. In 
other words, we should not expect the solution of (13. ip to belong to any (5)_p if p < 1. 

Theorem 3.2. Let = 0(x), x G M 6e a smooth compactly supported function and let 
^ be a standard Gaussian random variable. Then equation 

(3.4) Ut + uoux = Uxxi < t < T, a; G M, 

with initial condition u{0,x) = C,4'{x) cannot have a solution that is an element of 
(5)_p,_^(L2((0, T); H^{R)) for some 0<p<l andi>0. 

Proof. We will show that the ^-transform u oi u cannot be an entire function, as 
required by Theorem 12.61 for u to be an element of (iS)_p^, p < 1. 

By mi, 

(3.5) Ut(t, x; z) + u{t, x; z)ux(t, x; z) = Uxx{t, x; z), 

m(0, x; z) = z(j){x):, with only one random variable ^, we have only one complex param- 
eter z. 

Equation (13. 5p is the usual Burgers equation and has a closed-form solution via the 
Hopf-Cole transformation (see Evans [191 Section 4.4.1]): writing F{x) = J_ (j){y)dy, 
we find 

^, , rr#(y)exp(-(^-^)d, 

(3.6) u{t,x;z)= ^ ^ 



/r-p(-^-^)cii/ 



This is a classical solution of (13. 5p and leads to a classical chaos solution of (13. 4p : the 
solution of (13. 5p is unique in the class L2((0,T); H'^(M.)) (see Biler et al. [201 Theorem 
2.1]). If u{t, x) G (5)_p,_£(M) for every t, x, then, by Theorem 12.6( a). u{t, x; z) must be 
an analytic function of z for all t G [0, T], x G M, -2 G C. Since it is not immediately clear 
from (13. 6 p whether the dependence of m on 2; is analytic, we will transform equation 
(13. 6 p further by going back to the derivation of the Hopf-Cole transformation. 
Consider the equation 

(3.7) vt + -\vx\^ = Vxx, < t < T, xeR, 

with initial condition v{0,x;z) = zF{x) = z J^^(f){y)dy. Information about v leads 
to information about u because of the relation Vx = u. Note that F is a smooth 
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bounded function. Direct computations (see [191 Section 4.4.1]) show that the function 
V defined by 



(3.8) V{t,x;z)=expi-v{t,x;z)/2 
satisfies the heat equation 

(3.9) Vt = V;,, V^(0, x; z) = exp ( - zF{x)/2^ 

Equation (13. 9p has a classical solution which is unique in the class of bounded twice 
continuously differentiable functions. Interpreting the heat kernel as the normal den- 
sity, the solution of (13. 9p can be written as 

(3.10) \/(t,x;z)=Ee-"'^(*'"), 

where (,{^1^) = (1/2)-^ (a; + v^u;(t)) and tt; is a standard Brownian motion. 

We are now ready to show that the function u from (13. 6 p cannot be entire (that is, 
analytic for all z G C). The argument goes as follows: 

(1) The function V defined in (I3.10p is an extension to the complex plain of the 
characteristic function of a uniformly bounded non-degenerate random variable 
C, and is therefore an entire function of the form 

(3.11) V{t, x; z) = Vo{t, x; z)e^(*'^)^ 

where Vo{t,x; z) is an entire function with infinitely many zeroes (Lukacs [2T| 
Theorem 7.2.3]). 

(2) Representation (13. 8p then implies that the function v cannot be an entire func- 
tion of z; otherwise V would have no zeroes; 

(3) By uniqueness of solutions of (13. 5p and (13. 7p we conclude that Vxit,x;z) = 
u{t, x; z), because both functions satisfy the same equation with the same initial 
condition. 

(4) By the fundamental theorem of calculus, v{t, x; z) = v{t, 0; z) + J^ u{t, y; z)dy, 
which means that if u were an entire function of z, so would be v{t,x;z) — 
v{t,0;z). 

(5) If v{t, x; z) — v(t, 0; z) is entire, then, by (13.81) . so is 

t7/, n V{t,x;z) 

and moreover, V, as a function of z, has no zeros (because 1/V corresponds 
to v(t, 0; z) — v(t, x; z) and must be entire as well); then (13.111) and some basic 
complex analysis imply 

(3.12) V{t, x; z) = hit, x)e^^(*'^) or V{t, x; z) = V{t, 0, z)h{t, x)e"^(*'^\ 

(6) Finally, fl3.12p and the equality u = —2Vx/V imply that (if we assume u is 
entire, forcing v{t,x;z) — v{t,0;z) to be entire, in turn forcing V to have the 
form (I3.12P ) the solution u of (13. 5p must be a linear function of z, which is 
impossible and provides the required contradiction. 

Note that the key step in the argument, namely, showing that the function V(t, x; z) 
cannot have the same zeros in z for all x, can be carried out in several different ways. 
This concludes the proof of Theorem \3.^ D 
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Note that equation (13 ■4p without the Wick product, 
(3.13) t/i + f/f/, = f/,,, f/(O,x)=e0(a;), 

has a classical solution which is a very regular stochastic process. Indeed, the solution 
is given by the right-hand side of f l3.6p with ^ instead of z, from which it immediately 
follows that U{t, x) is a smooth function of (t, x) and 

for some positive numbers ci,C2. That is, as a random variable, U has moments of all 
orders. By contrast, the solution of (13. 4p (provided it exists) is a generalized random 
element from (iS)_i(M). In other words. Theorem (13. 2 p delivers the "curse" part of 
paper's title: Wick product in the Burgers equation indeed forces the solution into the 
largest space of generalized random elements. The following observation provides a bit 
of relief: with all the nice properties of U, there is no easy way to get KU{t,x); on 
the other hand, for the solution of fl3.4p . the generalized expectation u{t,x;0) solves 
the deterministic Burgers equation with zero initial condition and is therefore equal to 
zero. 

The next theorem delivers the "cure" part by establishing solvability of (13.1 p in 
(5)_i(L2((0,T); H'^{G))) For technical reasons, the conditions are slightly different for 
G = M and G = S^; in particular, we have to consider the "homogeneous" case / = 
when G = S^. We discuss these and other questions in the following section. Recall 
that, for a generalized random element $, E<l> denotes its generalized expectation, and 
E$ = $(0) = $(0). 

Theorem 3.3. Assume that 

• condition (13. 2 p holds; 

• Mo e {S)^i^^i{H\G)) for some i > and Euq G H\G); 

• If G = M, then also Emq G -^i(IR) and f(t,x) = Fx{t,x) for a function F such 
that both F and F^ are bounded and Holder continuous in {t,x); 

• If G = S^, then also / = 0. 

Under these assumptions, equation (13. ip has a unique solution and there exists an 
n > such that 

(3-14) MW-n:L2((0,T):HHR)) + SUp \\u{t, -^'t.^. mm < OO. 



-l,-n;L.((0,T);H^( „ ' ^^^^^^ '' 



It turns out that the generalized expectation M(o) of the solution has special signifi- 
cance. In particular, unlike the other chaos coefficients, M(o) and its partial derivative 
in X have to be bounded continuous functions, and the conditions of the theorem ensure 
that. This is yet another technical issue which we will discuss later. 

Example 3.4. Let us apply Theorem 13.31 to the Burgers equation with random vis- 
cosity 

(3.15) Ut + uou^= ([j2Q + W{t,x))ur,) , 

which was mentioned in the introduction as one of the motivations for investigating 
(13. ip . To concentrate on the effects of the noise in the viscosity, we assume that the 
initial condition is a non-random smooth function with compact support. 
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For the space-time white noise W{t,x), we have representation 
(3.16) Wit,x) = J2hkit,x)^k, 

k>l 

where {hk, A; > 1} is an orthonormal basis in L2((0,T) x G) (see, for example, Holden 
et ah pri, Definition 2.3.9]). Then fl3.15p becomes a particular case of ( 13. ip with 
Ofc = /io + hk, h = dhk/dx, Ck = gk = f = 0- 
If G = M, then we take 

hk{t,x) = hk,,k2{'t,x) = mk^{t)wk2{x), ki,k2>l, 

where mi are normalized sines and cosines, and wi are Hermite functions (the standard 
reference for Hermite functions is Hille and Phillips |22| Section 21.3]). In particular, 
it is known that sup^ |'W^/(a^)| < c/~^/^^ and therefore sup^. |w;(a;)| < d^^^"^ (recall that 
-w'{ + (1 + \x\'^)wi = 2lwi). Then (Q holds with qi = 21, r = 1. 
If G = S*^, then we take 

hk{t,x) = hki,k2{t,x) =mk^{t)mk2{x), ki,k2>l, 

where mi are normalized sines and cosines. In this case sup^. |?7iJ(x)| < cl, and again 
(|32D holds with qi = 21, r = 1. 

As a result, we can apply Theorem 13.31 with qi = 21 in both cases, which is the 
standard choice in the traditional white noise setting. 

We can generalize the model further and introduce, in addition to the white noise 
viscosity, a white noise random forcing and a white-in-space initial condition, with 
arbitrary correlation between the three (by suitably parsing the sequence $, for the 
construction of the corresponding random term). 

The eigenvalues of the type qi = (2/)^, r > 1 work in most models. The option to 
select other qi can help, for example, to deal with exotic stochastic forcing terms such 

as Zlfc {e^^sinx)^k- 

This concludes Example \3.4\ 

4. Outline of the proof and further directions 
Recall that our objective is to study the equation 

ut{t, x) + u(t, x) o u^(t, x) = u^^it, x) + fit, x) 



(^41) + ^(afc(t, x)u^^ + hk{t, x)ux + Cfe(t, x)u + gkit, x)) o ^k, 

k>l 

<t<T, X eG. 

By definition, the solution is a generalized random element with S'-transform u satis- 
fying 



(4.2) 



u{t,x; z) + / u{s,x; z)ux{s,x; z)ds = u{0,x; z) + / Uxxis,x; z)ds 

1.^1 "^ 
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The most natural approach to proving existence and uniqueness of solution for (14. ip 
is to prove existence and uniqueness of solution for fl4.2p . This approach is used, for 
example, by Benth et al. [I2] . Here are some of the reasons why we choose a different 
approach: 

• The complex- valued Burgers equation (14 ■2p is not easily handled by the Hilbert- 
space methods (a major reason is that, for complex numbers, z^ ^ kP), whereas 
our goal is to get solvability in the Sobolev spaces H'^ . 

• The Holder spaces could be a natural choice for (14. 2p . but the available methods 
lead either to existence that is local in time or require additional assumptions 
about the initial conditions (in the form of smallness of certain norms). We 
would like to avoid either of these restrictions. 

• The problem of uniqueness of solution for (14. 2 p is non-trivial. 

• Even if we succeed in finding u and proving that it is unique, there is still no 
easy way to get any regularity information about the solution, such as estimate 

Our approach is to get information about the solution from the chaos expansion. In 
what follows, we will use the notation 

(4.3) H«(0 = v^^c 
and an equivalent chaos expansion 

(4.4) u=Y,^cMc.{£.). 

By (I2.15p . representation (14.41) is more convenient for the computation of the Wick 
product. 

Our first step is to derive equations for the coefficients Ua- As a bonus, we also get 
an alternative way to prove uniqueness of solution. 

By (12. lip (and keeping track of the factorial terms), we get the following relation 
between the chaos coefficients Uq, from (14. 4 p and the S-transform u: 

(4.5) u{t,x,z) = y^^Ua{t,x)z°'. 

Before we proceed, let us introduce several notations. To reduce the number of sub- 
scripts, it is convenient to re-name the initial condition: Uq{x) = (p{x), and to use 
alternative notations for derivatives: ii = Ut, and Du = u^- With the generalized 
expectation M(o) having special significance (which also means frequent appearance in 
complicated formulas), we re-name it as well: M(o) = u. Now we are ready to derive 
the equations for Mq,. 
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Theorem 4.1. The function u is a (strong) solution of (14. 2p if and only if the chaos 
coefficients {u^, a G J^} satisfy (also in the strong sense) the system of equations 

(4.6) u + uDu = D^u, u(0,a;) = v9(o)(x); 

(4.7) Ue(fc) = D{Du^(k) - UMe(fc)) + a^D^ + hkDu + CkU + Qk, 
iie(fc)(0,a;) = v?e(fe)(x); 

(4.8) iio, = Di^Dua — uUq) — 2, UfiDuct-fB 

(0)<f3<a 

+ 2_^ [O'kD Ua-e(k) + bkDUa-e(k) + CkUa-e(k)) , 
fc:Ofc>0 

Ma(0,x) = ipa{x), \a\ > 1, 
and there exists a p > such that, for all ol ^ J , 

Proof Note that dS]) is equivalent to the condition u G (5)_i(L2((0,T); i/^^G))). 

Assume that u satisfies (14.21) and is an analytic function of z so that the series (14.50 
converges uniformly in {t, x, z) in some neighborhood of 2 = 0. We then substitute 
(14. 5 p into (14.21) . compare the coefficients of z° for each a, and get (I4.6l) - (l4.8p . Uniform 
convergence ensures that all manipulations are legitimate. 

For the proof in the opposite direction, we reverse the above argument. If (14.91) 
holds, then the series (14.51) converges uniformly in (t, x, z) in some neighborhood of 
z = 0. We then combine equalities (I4.6p - (l4.8p into power series and get (14. 2p . 
This concludes the proof of Theorem \4.1\ □ 



The system of PDEs (I4.6p - (l4.8p is called the propagator for equation (14. ip and 
describes how the stochastic equation propagates chaos through difference stochastic 
scales. Note that this propagation of chaos has no connection with the similar term 
used for deterministic equations in the study of particle systems (e.g. Sznitman [B]). 

Before we proceed, let us note that: 

(1) Since liuik^ooQk = +C)0, an upper bound q^" is equivalent to C'^'q^": in (14. 9p . 
we can always take a larger p, or, in the original construction of the chaos space, 
we can switch from qk to Cqk- 

(2) The propagator is a lower triangular system and can be solved by induction 
on \ck.\. Only the first equation, describing u = M(o) is non-linear: it is a 
deterministic Burgers equation. All other equations are linear, but with variable 
coefficients in the lower-order derivatives; these coefficients are determined by 
u. 

(3) By Theorem 14.11 uniqueness of solution for the propagator implies uniqueness 
of solution of (14. ip in the sense of Definition 13.11 

To prove Theorem l3.3l it is enough to show that the propagator has a unique solution 
and 

(4.10) lko|lL((o,T);//2(G)) + sup \\u„it, Wm{G) < c^'^'q''", p>0; 

0<t<T 

this is actually stronger than (13.141) . 
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A key step in the analysis of the propagator is the study of the Burgers equation 
(I4.6p . The solution of (14. 6p must be 

• a continuously-differentiable in x function, to act as a variable coefficient in the 
subsequent propagator equations, and also 

• an element of L2{{0,T); H'^{G)), to produce a sufficiently regular free term in 

62D. 

On the real line, the necessary solvability results can be obtained using the Hopf-Cole 
transformation, with sufficiently many regularity assumptions on ip and /. On the 
circle, we have to go with the best result we could find: Kiselev et al. [251 Theorem 
1.1], who consider / = 0. 

We will now present the precise results about solvability of the deterministic Burgers 
equation. 

Consider the equation 

(4.11) Utit, x) + f/(t, x)U.,it, x) = f/,,(t, x) + F,{t, x) 

with the initial condition f/(0, x) = Uo{x), for x E G where G = R (real line) or G = S*^ 
(unit circle or periodic boundary conditions). To define the solution of the equation, 
we assume that Uq is a continuous function and F G L2((0,T) x G)). 

A weak solution of (14. lip is a continuous in t, x function U such that, for all smooth 
functions tp with compact support in G, the equality 

U{t,x)il){x)dx — - I I U'^{s,x)il)'{x)dxdt = I UQ{x)il){x)dx 



(4.12) ''' ^ f '"^ /^ 

+ U{s,x)ip"{x)dxdt — / / F{s,x)'ijj'{x)dxds 

Jo Jg Jo Jg 

holds for allte [0,T]. 

Theorem 4.2 (Solvability of the deterministic Burgers equation). 

(a) If Uo G C^(R) f]Li(]R), Uq is Holder continuous (any non-zero order will work), 
and F, F^, are bounded and Holder continuous in (t, x) functions, then U-ll^ has 



a weak solution U = U{t, x) such that U is bounded and has a continuous bounded 
derivative in x. The solution is unique in the class of bounded continuous functions. 

(b) If Uq e C"^(M) f]-fr-'^(M) f]Li(R), U'q is Holder continuous (any non-zero order 
will work), and F and F^ are bounded and Holder continuous in {t, x) functions, then 
(^777]) has a weak solution U = U(t, x) such that 

(4.13) U e L,{{0,Ty,H\R))f^C{{0,Ty,H\R)). 

The solution is unique in the class of bounded continuous functions. 

(c) If Uq G H^{S^) and F = 0, then ^.11\ ) has a weak solution U = U{t,x) such 



that U G C((0, T); if^(S'"^)) (by the Sobolev embedding, this implies continuous differ- 
entiability in x). The solution is unique in the class of bounded continuous function. 

Proof. For the case G = M, existence, uniqueness, and regularity of the solution are 
established using the Hopf-Cole transformation: we have 

(4.14) U = -2VJV, 
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where V is the solution of the heat equation 

(4.15) Vt = V,,-^F{t,x)V, r(0,x) = exp(-^y" f/o(y)rfi/). 

Thus, the study of equation (14.121) is reduced to the study of equation (I4.15p . 

(a) A result from Ladyzhensakaya et al. p3| Theorem IV. 5.1] implies that, under 
our assumptions, equation (I4.15P has a classical solution and the solution is unique in 
the class of bounded continuous functions. Then a result from Rozovskii |2H Theorem 
5.1.1] provides the following probabilistic representation of the solution, with w = w{t) 
denoting a standard Brownian motion: 

(4.16) V{t,x)=E(Vo{x + V2w{t))exp{-- / F{s,x + V2w{t) - V2w{s))ds) 

Since F is bounded and 

g-(i/2)||c/ollLi(R) < V(x) < e(^/2)ll^oll^iW, 

it follows from (I4.16P that there exist numbers < pi < p2 such that, for all t, x, 
Pi < V{t, x) < p2. Therefore the functions U = -2V^/V and U^ = -2{V^^V - V^)/V^ 
are bounded and continuous. Uniqueness of U follows from the uniqueness of V. The 
reader can also note that considering the Burgers equation in the form Ut+UUx = Uxx/2 
would eliminate a lot of various powers of 2 in subsequent computations. 

(b) Under conditions of this part, equation (I4.15P still has a classical solution, but 
the standard parabolic regularity theorem in Sobolev spaces cannot be applied to V 
because Vq is not integrable, and so V does not belong to any Sobolev space. On the 
other hand, (I4.14p implies that (I4.13P will follow if we can show that 1/V is bounded 
and has two continuous and bounded derivatives in x (which we did in part (a)), and 
that 

(4.17) \4GL2((0,T);i72(R))f|C((0,T);ifi(M)), 

which we show next. 

Denoting Vx by V, we find from (I4.15P that 

Vt = Vxx-^FV-^VFx, 

and V{0, x) = —{l/2)Uo{x)Vo{x). Assumptions of the theorem and the properties of V 
imply that V^(0, ■) G H^ and VFx € -^2((0, T) x R). Then the linear parabolic regularity 
theorem (see Theorem 14.31 below) in the normal triple {H'^{M.),H^(M.),L2{M.)) implies 

mm. 

By the Sobolev embedding, the conditions Uq G C^(M) f] H^{M.) plus Holder continu- 
ity of Uq can be replaced with a stronger but shorter condition Uq G H'^{M.), and this 
is what happens in the statement of Theorem 13.31 In general, we do have to add the 
condition Uq G -Li(M) (for example, the function (1 + |xp)~^'^ is in every if'*'(M), but 
is not in Li(M).) Finally, note that the resulting weak solution U, being an element of 
L2{{0,T);H'^{R)), is in fact strong. 
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(c) See Kiselev et al. [251 Theorem 1.1]. 
This completes the proof of Theorem \4-^ D 

We were unable to find satisfactory solvability results for the circle when an inho- 
mogeneous term f{t,x) is present, or for the Dirichlet or Neumann boundary value 
problems. As soon as these results are available, the corresponding analogue of The- 
orem [33] will follow immediately from Theorem 14.11 In fact, one can go even further 
and investigate Wick-stochastic versions of other one-dimensional equations with qua- 
dratic nonlinearity, such as Camassa-Holm, KdV, KPP, and Kuramoto-Sivashinskii 
equations, provided there is a satisfactory solvability result for the corresponding deter- 
ministic nonlinear equation. 

Extension to higher dimensions presents an additional challenge, as the Sobolev space 
H^(R'^), d > 2, is no longer embedded into the space of continuous functions. Never- 
theless, Mikulevicius and Rozovskii [TB] recently studied the Wick-stochastic versions 
of the Navier-Stokes equations using Lp solvability results with p > d. Extension of 
our results to a (i- dimensional system of Burgers equations seems feasible, especially on 
W^, where the Hopf-Cole transformation reduces the system to a scalar heat equation. 

Similar to [16], one can study other properties of the chaos solution of (particular 
cases of) (14.11) . such as adaptedness and Markov property. For a number of reasons, 
these questions fell outside the scope of the current paper (in fact, to simplify presen- 
tation, we do not even consider a filtration and introduce random perturbation simply 
as a countable set of iid standard Gaussian random variables). 

We now continue with an outline of the proof of Theorem l3.3[ The difference between 
G = R and G = S^ stops with Theorem 14. 2[ After that, analysis of the propagator is 
all about linear parabolic equations of the form 

ut = Au + F, 

where ^ is a linear partial differential operator. The analysis relies on the following 
result. 

Theorem 4.3 (Linear Parabolic Regularity). Let {V,H,V') be a normal triple of sep- 
arable Hilbert spaces and denote by [w,v], w G V, v eV, the duality between V and 
V relative to the inner product in H . Consider a collection of bounded linear operators 
•A{t), <t <T, from V to V with the following properties: 

• For all x,y E V , the function t i-> [A{t)x, y] is measurable; 

• There exist e > and Cq G M such that, for all t G [0, T] and v E V, 



[A{t)v,v]+6\\v\\l < Co\\v\ 



H- 



Then, for every uq E H and F G L2((0,T); V), there exists a unique u G L2((0,T); V) 
such that the equality 

u{t) = UQ+ I A{s)u{s)ds+ I F{s)ds 
Jo Jo 
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holds m L2{{0,Ty,V'). Moreover, u G C((0,T);i7) and 

[ \\u{t)\\ldt+ sup \\u{t)\\l < C{e,Co,T){\\uo\\l + f \\F{s)fy,ds). 

Jo 0<t<T Jo 

Proof. See, for example, Rozovskii pH Theorem 3.1.2]. D 

To study the propagator equations (14.81) . we apply Theorem 14.31 with 

• V = H\G), H = H\G), V = L2(G); 

• A = d'^/dx'^ — u(t, x)d/dx — Ux{t, x), where u is the solution of (14.61) : 

• A complicated but manageable free term 

F = Fc = - ^ Ui3DUa^i3 + ^ {akD^Uc-e(k) + bkF>Uc-e(k) 
(4.18) (0)<f3<a k:ak>0 

+CkUa-e(k)) ■ 

By Theorem 14. 2[ the functions u and Ux are continuous and bounded, so that the oper- 
ator A is indeed bounded from H'^{G) to L2{G). The terms Ujs have (3 < a. and come 
from the earlier propagator equations. In particular, by Theorem 14.31 and the Sobolev 
embedding, Ufj G L^((0,T); if^(G)) f]C((0, T) x G), which means that expressions 
Uf^Dua-fs and akD'^Ua-e(k) are square-integrable, and therefore F G L2((0,T) x G). 

The main technical complication is that the number of terms on the right-hand side 
of (I4.18P grows with \a.\, and, to satisfy (l4.1Up . we have to control this growth. In the 
following section, we derive a recursive relation between the L2 norms of F^ for different 
ex. and show that the rate of growth in \a.\ is controlled by the Catalan numbers C\ct\-i 
(along with the admissible factor q^"). The exponential growth rate of the Catalan 
numbers {Ck < 4^^) allows us to complete the proof. 

5. Analysis of the propagator 

In this section we carry out the analysis of the propagator (I4.6p - (I4.8I) and complete 
the proof of Theorem 13. 3[ Our main goal is to establish (I4.10p . 

Equation (14. 6 p is covered by Theorem 14.21 under the assumptions on y9(o) = Emq 
and /, we know that u and u^ are bounded and continuous, and 

ueL2iiO,Ty,H\G)). 

Next, consider (14.70 . If ck = e^, we have 

«6fe =D{DUe, -M(o)MeJ + F^, (t, x), U^, (0, x) = ^e^, 

where 

Fe^(t,x) = afcD^M(o) + bkDu(o) + CkU(o) + 9k 

By Theorem 14.21 F^^ G L2 ((0,T) x G) . We then solve for u^^. in the normal triple 
{H'^{G),H\G),L2{G)) using Theorem SJ and deduce fHTTUD for |a| = 1. 

As was mentioned earlier, induction on \a\, Sobolev embedding, and Theorem 14.31 
imply that F^, the free term in (14. 8 p (see also (I4.18P ). is an element of L2((0,T) x G). 
All we need is a bound on the norm of Fa that is consistent with (I4.10p . 



20 S. KALIGOTLA AND S. V. LOTOTSKY 

By the Sobolev embedding, the L2 norm of F^ with \cx\ > 1 is controlled by the norms 
oiu,3, (0) < /3 < a, in the space L2{{0,T); H^{G)) f]C{{0,T); H\G)). Accordingly, 
for ex. with \a\ > 1, we define 

-^a = / \\'^cx{tr)\\2dt + sup \\Ua{tr)\\l = \\Ua\\L2((0.T)-,m(G)) ~^ \\'^ci\\c({0,Ty,m(G))- 
Jo 0<t<T 

The next step is to derive a recursion for L^- 
By Theorem 14.31 applied to (14. 8p , 

(5.1) ^o<A(||</^o||?+||F«||i,((0,T)xG))> 

and A (here and below) is a number depending only on T and some norms of u. The 
value of A can change from line to line. With no loss of generality, we assume that 
A > 1, and, with m = max(£, r), put 

Le(k) = MT- 

For \(y.\ > 1, (15.11) . fl4.18p . and assumptions of Theorem 13.31 imply 

r ^ \ ( f^ct \ ^ II II II II 

-t^a < Alq + 2_^ \\Ul3\\c{[0,T]xG) ||'"a-/3 ||L2((0,T);/fi(G)) 

(0)</3<a 

fc:afe>0 

< A 



V (0)</3<Q k:ak>0 / 



Define 



Then 



L, = 2A ( ^ + 1 1 . 



(0)<f3<a 

Let {Aa : a G JT"} be a sequence such that A^^, = I/e(fc) and 

(5.2) Aa= Yl ^/3^o-/3, |a| > 1. 

(0)</3<a 

It follows by induction that 

(5.3) La < Aa 



for all OL ^ J . Indeed, (15. 3p is true by definition if \ol\ = 1. Assume (15. 3p holds for all 

|a| < A^. For \cx\ = N, 

La < 2_^ L/^La-lB 

(0)</3<a 

< Y, ^I3^»-I3 = Aa, 

(0)</3<a 
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proving (15.31) for all a G JT. 

Our next step is to find a bound on Aa, which we do using generating functions. 
Since multi-indices o: can have non-zero entries in arbitrary positions, the argument 
requires an additional construction. 

Let d > 1 he an integer and let /3 = {f3i, (32, . . . , f3d) be a c?— dimensional vectors 
of nonnegative integers with at least one /3j > 0. Let M^. G M, i = 1,2, ... ,d, be 
prescribed positive numbers and let Mf3, |/3| > 1, be defined through the recursive 
equation, 



7+<5=/3,|7|,|(5|>0 



Consider, 



F,{z^,Z2,....,z,) = J2K4 



id) ^/3i 82 e 



Then 



Z2 ■ ■ -z^ 



FJ{z„Z2,....,z,) = Y.\ E <^Mf hf^ 

|/3|>1 \7+5=/3,|7|,|5|>0 / 

i/3i>i 

i 

where e^ is the (i- dimensional vector with 1 at the i*'* coordinate and everywhere else. 
Solving for F^ we get. 



Fd[Zi, Z2, ...., Zd) = 

By comparing the coefficients of similar powers we obtain. 



(5.4) 



<'^M(r:0('^')n«''' 



here (';^') = \/3\l//3\. Let L^ = {ot = {ai,a2, . ■ ■) e J : ai = for all i>d}. If 
a E Td then, by (15. 2p . A^ is uniquely determined by {Af^ : j3 < a.}. Also, if a G F^ 
and (3 < oc, then /3 G Ff^ as well. Set MeJ = A,^., for each e^ G F^. Then 



(5.5) M^l = A 



a. 



where cx'^ = (ai, 0^2, ... , ad), which follows by induction on \a\. Indeed, equality (15.51) 
is true if a G F^ such that |a| = 1. Assume (15.51) is true for a G F^ such that \a\ < N. 
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For a G r^ such that \a\ = N, 

7+5=«,|7|,|<5|>0 

5] Tl^l^M 



id) T^/fid) 
7d+5'*=a'*,|7'*|,|5'*|>0 



and the general equahty fIS.Sp follows. 

Given any ol & J there exists a (i G N such that a G F^. For such a (i, M^ j = Aa,. 
Using ([52D, 

. N . 1 /2|a| -2\ f\a.\\ tt ,„ 

q; V a - 1 / V a /-'■-'■ ^ 



We now notice that 



1 /2r2\ 

Cn = —— 1 , n > 0, 

n + 1 \ n 



is nth Catalan number, a very popular combinatorial object (see, for example, Stanley 
[26]). Hence, re-tracing our arguments back to L^, we conclude that 



(5.7) 



i:c.<(2A)2Hq^i_/H'jq^ 



The final two observations are 

• C*n < 4", which follows by the Stirling formula; 

• (") < l") which we proved in Proposition 12.31 

Therefore, L^ < (8A)'"'q^™'^^^", which establishes (I4.10p and completes the proof of 
Theorem 13.31 
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